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Abstract 

Nonlinear charge transport in strongly coupled semiconductor superlattices is described by 
Wigner-Poisson kinetic equations involving one or two minibands. Electron-electron collisions 
are treated within the Hartree approximation whereas other inelastic collisions are described by a 
modified BGK (Bhatnaghar-Gross-Krook) model. The hyperbolic limit is such that the collision 
frequencies are of the same order as the Bloch frequencies due to the electric field and the corre- 
sponding terms in the kinetic equation are dominant. In this limit, spatially nonlocal drift-diffusion 
balance equations for the miniband populations and the electric field are derived by means of the 
Chapman-Enskog perturbation technique. For a lateral superlattice with spin-orbit interaction, 
electrons with spin up or down have different energies and their corresponding drift-diffusion equa- 
tions can be used to calculate spin-polarized currents and electron spin polarization. Numerical 
solutions show stable self-sustained oscillations of the current and the spin polarization through a 
voltage biased lateral superlattice thereby providing an example of superlattice spin oscillator. 

PACS numbers: 



I. INTRODUCTION 

Semiconductor superlattices are essential ingredients in fast nanoscale oscillators, quan- 
tum cascade lasers and infrared detectors. Quantum cascade lasers are used to monitor 
environmental pollution in gas emissions, to analyze breath in hospitals and in many other 
industrial applications^. A superlattice (SL) is a convenient approximation to a quasi- 
one-dimensional crystal that was originally proposed by Esaki and Tsu to observe Bloch 
oscillations, i.e., the periodic coherent motion of electrons in a miniband in the presence of 
an applied electric field. Fig. [It^a) shows a simple realization of a A^-period SL. Each period 
of length / consists of two layers of semiconductors with different energy gaps but with sim- 
ilar lattice constants. The SL lengths in the lateral directions, Ly and L^ are much larger 
than /, typically tens of microns compared to about ten nanometers. The energy profile of 
the conduction band of this SL can be modeled as a succession of square quantum wells 
and barriers along the x direction (Kronig-Penney model) and, for a n-doped SL, we do not 
have to consider the valence band. A different quasi- ID crystal called a lateral superlattice 
(LSL) is shown in Fig. [U^b). In this case, a periodic structure is formed on the top surface 
of a quantum well (QW), so that L^ is of the order of / and Ly ^ /. The wave functions 
of a single electron in the conduction band of a SL can be expanded in terms of ID Bloch 
wave functions times plane waves 

^ '^yyi){z)e'^''u^{x,k), (1.1) 

, for a SL, 

(1.2) 
tpniz), for a LSL, 

where u is the miniband index and n is the energy level of the quantum well in the case of 
a LSL. The function u^{x, k) is /-periodic in x and 27r//-periodic in k. S is the area of the 
lateral cross section, equal to LyL^ for a rectangular cross section. 

Many interesting nonlinear phenomena have been observed in voltage biased SL com- 
prising finitely many periods, including self-oscillations of the current through the SL due 
to motion of electric field pulses, multistability of stationary charge and field profiles, and 
so on^. It is important to distinguish between strongly and weakly coupled SLs depending 
on the coupling between their component QWs. Roughly speaking, if barriers are narrow, 
QWs are strongly coupled and we can use the electronic states (11.11) as a convenient basis in 
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FIG. 1: (a) Schematic drawing of a superlattice. (b) A lateral superlattice. 



a quantum kinetic description. The resulting reduced balance equations for electron density 
and electric field are partial differential equations (which may be nonlocal, as we shall see in 
this paper). On the other hand, for SLs having wide barriers, their QWs are weakly coupled 
and the electronic states of a single well provide a good basis in a quantum kinetic descrip- 
tion, replacing the Bloch functions e^^^u^[x^ k) in (11. ip . In this case, the balance equations 
are spatially discrete and phenomena such as multistability of stationary field profiles, for- 
mation and pinning of electric field domains, etc are theoretically predicted and observed in 
experiments. See the review^. Another promising field of applications is spintronics. Elec- 
trons in SLs having at least one period doped with magnetic impurities and subject to a 
static magnetic field can be distinguished by their spin because the magnetic field splits each 
miniband in two having different spin-dependent energy—. Recently a SL of this type has 
been proposed as a spin oscillator producing spin polarized oscillatory currents and able to 
inject polarized electrons in a contact^. Alternatively, materials displaying strong spin-orbit 
effects can be used as spintronic devices without having to apply magnetic fields; cf. the 
case of the LSL considered in Ref.— . In this paper, we will show that a LSL can be used as 
a spin oscillator. 

This paper presents systematic derivations of quantum balance equations for SLs with 
two populated minibands, and it shows that their numerical solutions may predict space 
and time-dependent nonlinear phenomena occurring in these materials. Our methods can 
be used in 3D crystals, but their application to ID structures such as SLs and LSLs leads to 
simpler equations that are less costly to solve. Although nonlinear charge transport in SLs 
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has been widely studied in the last decade (see the reviews^'^^'^i) ^ systematic derivations 
of tractable balance equations for miniband populations and electric field are scarce. One 
reason is that quantum kinetic equations are nonlocal in space and their collision terms 
may be nonlocal in space and time^'^. Using them to analyze space and time- dependent 
phenomena such as wave propagation or self-sustained oscillations is problematic. In fact, 
only extremely simple solutions of general quantum kinetic equations (such as thermal equi- 
librium, disturbances thereof due to weak external fields and so on) are known, theoretical 
analysis of these equations is lacking and numerical solutions describing spatio-temporal 
phenomena are not available. One way to proceed is to adopt simple collision models sim- 
ilar to the Bhatnagar-Gross-Krook (BGK) collision model for classical kinetic theory^. We 
discuss in this paper how to implement a BGK collision model for a quantum kinetic equa- 
tion that is simple to handle yet keeps an important quantum feature such as the broadening 
of energy levels^. Once we have a quantum kinetic equation for a sufficiently general SL 
having two minibands, we implement a Chapman- Enskog perturbation procedure to derive 
the sought balance equations and solve them numerically for realistic SL configurations. 

Previous to this work. Lei and coworkers derived quantum hydrodynamic equations de- 
scribing SL having only one miniband^»^. They use a closure assumption to close a hi- 
erarchy of moment equations. For the case of quantum particles in an arbitrary external 
three-dimensional potential, Degond and Ringhofer- have used a similar procedure to derive 
balance equations. They close the system of moment equations by means of a local equilib- 
rium density obtained by maximizing entropy. The Chapman-Enskog method has been used 
to derive drift-diffusion equations for single-miniband SLs described by semiclassical^ and 
quantum kinetic equation^. Earlier, Cercignani, Gamba and Levermore used the Chapman- 
Enskog method to derive balance equations for a semiclassical BGK-Poisson kinetic descrip- 
tion of a semiconductor with one parabolic band under strong external bias^. 

The rest of this paper is as follows. In Section [TTl we review the simpler case of nonlinear 
electron transport in a strongly coupled n-doped SL having only one populated miniband^. 
Starting with a kinetic equation for the Wigner function, we use the Chapman-Enskog per- 
turbation method to derive balance equations for the electron density and the electric field. 
When these equations are solved numerically for a dc voltage biased SL with finitely many 
QWs and realistic parameter values, stable self-sustained oscillations of the current through 
the SL are found among their solutions, in agreement with experimental observational. Sec- 
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tions mil to |V] contain the main results of the present work. In Section llllt we describe a SL 
having two populated minibands by proposing a kinetic equation for the Wigner matrix. In 
Section IIVI we derive balance equations for the miniband electron populations and the elec- 
tric field, using an appropriate Chapman-Enskog method and a tight-binding approximation 
to obtain explicit formulas. The case of a LSL having strong Rashba spin-orbit interaction^^ 
is important for spintronic applications and has been considered in Section |Vl We derive 
and solve numerically the resulting balance equations. Novel self-sustained oscillations of 
the spin current and polarization are obtained for appropriate values of the parameters. 
Finally Section IVII contains our conclusions and some technical matters are relegated to the 
Appendix. 

II. SINGLE MINIBAND SUPERLATTICE 



The Wigner-Poisson-Bhatnagar-Gross-Krook (WPBGK) system for ID electron trans 
port in the lowest miniband of a strongly coupled SL is: 
df . t 



dt h 

ie 
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(2.4) 



knT J\ [E-£^ik)]^ + T^ 
Here /, n, N^, £ik), ds, dw, I = ds + dw, W, e, m*, ks, T, T, u^n, t'imp and — e < 
are the one-particle Wigner function, the 2D electron density, the 2D doping density, the 
miniband dispersion relation, the barrier width, the well width, the SL period, the electric 
potential, the SL permittivity, the effective mass of the electron in the lateral directions, 
the Boltzmann constant, the lattice temperature, the energy broadening of the equilibrium 
distribution due to collisions'^ (page 28 ss), the frequency of the inelastic collisions respon- 
sible for energy relaxation, the frequency of the elastic impurity collisions and the electron 
charge, respectively. 
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The left-hand side of Eq. (12 .ip can be straightforwardly derived from the Schrodinger- 
Poisson equation for the wave function in the miniband using the definition of the ID Wigner 
function^: 



f{x,k,t) 



2/ ^ 



s 



j=-oo 



{^\x + jl/2, y, z, t)^{x - jl/2, y, z, t))e^^'^'rf: 



x± 



(2.5) 



(the second quantized wave function ^(x, x_L,t) = X^gq (^iQ^ Q±^'^)4'qi^)^^^^^^ ^ ^± = iu^z): 
is a superposition of the Bloch states corresponding to the miniband and S is the SL cross 
section-). The right hand side in Eq. (12. ip is the sum of — z/g (/ — f^^), which represents 
energy relaxation towards a ID effective Fermi-Dirac (FD) distribution f^^{k; n) (local equi- 
librium), and —Vi[f{x,k^t) — f{x, —k,t)]/2, which accounts for impurity elastic collisions^. 
For simplicity, the collision frequencies z/g and Ui are fixed constants. Exact and FD dis- 
tribution functions have the same electron density, thereby preserving charge continuity as 
in the classical BGK collision models^. The chemical potential /i is a function of n result- 
ing from solving equation (12. 3p with the integral of the collision-broadened 3D Fermi-Dirac 
distribution over the lateral components of the wave vector (/c, k_|_) = {k, ky, k^): 
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Using the residue theorem for a line-width: 
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which is equivalent to Eq. (I2.4p 2^. Here 0[E) is the Heaviside unit step function. As 
F -^ 0+, the hue- width ([23]) tends to the delta function b{E), Dr{E) tends to the 2D 
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density of states, D{E) = 'm*6{E)/{TTfi^), and f^^ tends to the 3D Fermi-Dirac distribution 
function integrated over the lateral wave vector k_L. In Ref.-, a Lorentzian line- width was 
used instead of fl2.8p and the integral over E in (12.61) extended from to cxo. The integral 
with the Lorentzian function is not convergent in i^^ = — oo, which is why we prefer using 
convolution with the "super-Lorentzian" function fl2.8p in this work. The integration in 
(12. 7p cannot be carried out explicitly for other standard line-widths such as a Gaussian 
or a hyperbolic secant. This unnecessarily complicates the numerical integration of the 
balance equations we will obtain later. Note that, following Ignatov and ShashkinW, we 
have not included the effects of the electric potential in our Fermi-Dirac distribution. These 
model equations can be improved by including scattering processes with change of lateral 
momentum and an electric field-dependent local equilibrium. However the resulting model 
could only be treated numerically and the qualitative features of our derivation and of the 
nonlocal drift-diffusion equation would be lost in longer formulas. 

A different way to introduce a quantum BGK collision model is to define a local equi- 
librium density matrix operator by minimizing quantum entropy (defined with the opposite 
sign of the convention that is usual in physics) under constraints giving the electron den- 
sity and energy density in terms of the density matrix. The resulting expression involves 
an inverse Wigner transform and another transform is needed to deduce the local equilib- 
rium Wigner function f^^ entering the BGK formula^. This f^^ is nonlocal in space and 
can only be found by solving some partial differential equation^. As a model for quantum 
coUisions^i^, the resulting quantum BGK model is not realistic, in the same way as the 
original BGK model is not a realistic model for classical collisions. Moreover, the implicit 
manner in which the model is defined defeats the main asset of the classical BGK collision 
model: its simplicity, that makes it possible to obtain results analytically. Thus we prefer 
to introduce a BGK model that can be handled more easily and still incorporates quantum 
effects. The most important quantum effect affecting the collision term is the broadening 
of energy levels due to scattering, F ^ h/r (where r is the lifetime of the level)^, and this 
is taken phenomenologically into account by the convolution with the line-width function 
(12. 8p in (12. 6p . In the semiclassical limit "/i — > 0", F — ;> and we recover the semiclassical 
FD distribution. 

The WPBGK system (12. ip to (12. 4p should be solved for a Wigner function which is 27r//- 
periodic in k and satisfies appropriate initial and boundary conditions. It is convenient to 
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derive the charge continuity equation and a nonlocal Ampere's law for the current density. 
The Wigner function / is periodic in k; its Fourier expansion is 

oo 

fix,k,t)= Y, fA^,t)e'^'K (2.10) 

j=-oo 

Defining F = dW/dx {minus the electric field) and the average 

it is possible to obtain the following equivalent form of the Wigner equation^ 



(F)Jx,t) = — F(x + s,t)ds, (2.11) 

jl J-il/2 



j=-oo 

Here the nonzero Fourier coefficients of the dispersion relation are simply Sq = A/2 and 
S±i = —A/4 for the tight-binding dispersion relation S{k) = A (1 — coskl)/2 (A is the 
miniband width), which yields a miniband group velocity v{k) = ^sin kl. Integrating this 
equation over k yields the charge continuity equation 

'i^r.tfi^'^i^-Mh-o- P.13) 

Here we can eliminate the electron density by using the Poisson equation and then integrate 
over X, thereby obtaining the nonlocal Ampere's law for the total current density J{t): 

O 77' O '^^^ 

i=i 
To derive the QDDE, we shall assume that the electric field contribution in Eq. (I2.12p is 
comparable to the collision terms and that they dominate the other terms {the hyperbolic 
limit)^. Let vm and Fm be the electron velocity and field positive values at which the (zeroth 
order) drift velocity reaches its maximum. In this limit, the time to it takes an electron with 
speed Vm to traverse a distance Xq = eFMl/{eND), over which the field variation is of order 
Fm, is much longer than the mean free time between collisions, z/^^ ~ h/{eFMl) = ti. We 
therefore define the small parameter X = ti/to = hvnfNu/{eFljl'^) and formally multiply the 
first two terms on the left side of (12.11) or (I2.12p by A-^^. The result is 



oo 



M w + E X ^'"^' hf)' = w\ - E ¥ '^""(P)' !'■ (2-15) 



dt ^^-^ h dx ^-^ h 

-oo / jf=— oo 



NONLINEAR ELECTRON AND SPIN TRANSPORT IN SUPERLATTICES 10 



The solution of Eq. (12.151) for A = is calculated in terms of its Fourier coefficients as 



f^'\k; F) = f; ^^ i'^tlV'"'' e--^ (2.16) 



j=-oo 

where J^j = {F)j/Fm, Fm = ^^^vj^v^ + z/j) and t^ = \/{ve + Vi)lve 
The Chapman-Enskog ansatz for the Wigner function is^: 



/(x, fc, t- A) = j^\k- F) + J2 f'^Kk- F) X^, (2.17) 



m=l 



e^ + Yj^'^\F)\'^ = J{t). (2.18) 

ot ^-^ 

The coefficients f^"^\k;F) depend on the 'slow variables' x and t only through their de- 
pendence on the electric field and the electron density. The electric field obeys a reduced 
evolution equation (I2.18P in which the functionals J^"^\F) are chosen so that the f^"^\k; F) 
are bounded and 2'7r//-periodic in k. After we keep the desired number of terms and set 
A = 1, Eq. ( 12.18P is the QDDE provided by our perturbation procedure. 

Differentiating the Ampere's law (12.181) with respect to x, we obtain the charge continuity 
equation. Moreover the compatibility condition 

f^"'\k-n)dk=^ft^ = 0, m>l, (2.19) 

■w/l '■ 

is obtained by inserting the expansion (I2.17P into (12.31) . Inserting (12.171) and (12.181) in (12.151) . 
we find the hierarchy: 






(1) _ _^ 



dt 



^ E ^^if-^^ (-0) 



j=— oo 



dt 
and so on. Here 



^ E ^lif'% - 1/'°' 



j=— oo 



(2.21) 



oo 

Cu{k) ^ '^Y.jl{F),u,e^'^ + (v, + I) u{k) - I u{-k), (2.22) 

— oo 

and the subscripts and 1 in the right hand side of these equations mean that e dF/dt is 
replaced by J — J^^\F) and by —J^^^F), respectively. 
The condition (I2.19P implies that 

/ Cf^'^Uk = 0, (2.23) 

J-n/l 
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for m > 1. Using this, the solvabihty conditions for the hnear hierarchy of equations yield 

^^'"^ = |E^'«^-^/r^))^' (2.24) 

i=i 

which can also be obtained by insertion of Eq. (12.171) in fl2.14p . 

Particularized to the case of the tight-binding dispersion relation and F = in the Fermi- 

Dirac distribution (12.41) . the leading order of the Ampere's law fl2.18p is 

^% + ^i^^y(^))i = -J(t), (2.25) 

/TV 
cos(mfc) In (1 + e^-5+5cosfc^ ^^^ ^2.27) 

■TT 

provided JF = JF^, S = A/(2A;bT) and fi = ^/{ksT). Here M (calculated graphically in 
Fig. 1 of Ref.-) is the value of the dimensionless chemical potential fi at which (12. 3p holds 
with n = No- The drift velocity vmVIJ-") has the Esaki-Tsu form with a peak velocity that 
becomes vm ~ A//i(5)/[4/ire/o((5)] in the Boltzmann limiti^ {In{S) is the modified Bessel 
function of the nth order). 

To find the first-order correction in (I2.18p . we first solve (12.200 and find J^"^^ for m = 1. 
The calculation yields the first correction to Eq. (I2.25P (here ' means differentiation with 
respect to n)^: 



1 



M = {nVM), + {{A - l){{nVM),),), - ^\ (t4^) ' (2-30) 

8n2(z/e + z/,)(l+-^') V dx^ Al J' ^' ' 
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^T^nM^ d{F)A ^/ nM,{l-A:Fl)d{F)A 
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AhvMil + rDJ'inM)' /^^ 1-^2 Q^p^^ 



A/re(l + J^2) \ (1+J^2)2 Q^ 



_ / {nM,y d'F\ _ /{uM^md'FX 
^ " \l + 4J-| dx^/, ^\ 1 + 4^1 dx^/, ^ "^"^^ 

8;it;M(l + T^){nMyj^ / {nUyT d^F 



A/re(l+J^2^ \ 1 + J^2 5^2/^ 
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Here M2{n/ND) = l2{fl) Tq{M) / [Xi{M) Xo^fl)]. If the electric field and the electron density 
do not change appreciably over two SL periods, {F)j ^ F, the spatial averages can be ig- 
nored, and the nonlocal QDDE (12 .28^ becomes the local generalized DDE (GDDE) obtained 
from the semiclassical theoryS'. The boundary conditions for the QDDE f l2.28p (which con- 
tains triple spatial averages) need to be specified on the intervals [—2/, 0] and [Nl, Nl + 21], 
not just at the points x = and x = Nl, as in the case of the parabolic GDDE. Similarly, the 
initial condition has to be defined on the extended interval [—21, Nl + 21]. For realistic values 
of the parameters representing a strongly coupled SL under dc voltage bias, the numerical 
solution of the QDDE yields a stable self- sustained oscillation of the current^ in quantitative 
agreement with experiments^. Details of the numerical procedure can be found in^. 

III. WIGNER DESCRIPTION OF A TWO-MINIBAND SUPERLATTICE 

We shall consider a 2 x 2 Hamiltonian H(x, —id/dx), in which^ 

H(x, k) = [ho{k) - eW{x)](7o + h{k) ■ a] 

{a + 'y){l — cos kl) — eW{x) + g —iPsinkl 

iPsinkl {a — 'y){l — cos kl) — eW{x) — g 



■ (3.1) 



Here 



ho{k) = a {1 — cos kl) , hi{k) = 0, 

h2{k) = P sin kl, hs{k) =7(1 — cos kl) + g. 



(3.2) 



and 





(T0= \ ' \,(T,= \ ' \,CT2= [^ M , 0-3 = K ° I (3.3) 

are the Pauli matrices. 

The Hamiltonian (13.11) corresponds to the simplest 2x2 Kane model in which the 
quadratic and linear terms {klY/2 and kl are replaced by (1 — cos kl) and sin kl, respec- 
tively. For a SL with two minibands, 2g is the miniband gap and a = (Ai + A2)/4 
and 7 = (Ai — A2)/4, provided Ai and A2 are the miniband widths. In the case of a 
LSL, (7 = 7 = 0, and h2(T2 corresponds to the precession term in the Rashba spin-orbit 
inter act ion?!^. The other term, the intersubband coupling, depends on the momentum in 
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the y direction and we have not included it here. Small modifications of (13. ip represent a 
single miniband SL with dilute magnetic impurities in the presence of a magnetic field B: 
(7 = 7 = /i2 = 0, and hi = [3{B)^. As in the case of a single miniband SL, W{x) is the 
electric potential. 

The energy minibands £^{k) are the eigenvalues of the free Hamiltonian Ho(fc) = 
hQ{k)(TQ + h{k) ■ <T and are given by 



S±(k) =ho{k)±\h{k)\. 



(3.4) 



The corresponding spectral projections are 



p±(^) ^ o-o±y)-o- ^ ^^^^^ ^^^^ ^ ^^^^^i^^^^i^ 



(3.5) 



so that we can write 



Ho(A;) = S+{k)P+{k) + S-{k)P~{k). 



(3.6) 



We shall now write the WPBGK equations for the Wigner matrix written in terms of the 
Pauli matrices: 

3 



f (x, k,t) = Y^ fix, k, t)(Ti = fix, k, t)(To + fix, k, t) 



(3.7) 



i=0 



The Wigner components are real and can be related to the coefficients of the Hermitian 
Wigner matrix by 



/21=/^+^/^ f22 = f-f. 



(3.8) 



Hereinafter we shall use the equivalent notations 



/ 



f 



p 
v/v 



(3.9) 



The populations of the minibands with energies £^ are given by the moments: 



n ix,t) 



27r 



■k/I 



fix,k,t)±Vik)-fix,k,t) 



dk. 



(3.10) 
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and the total electron density is n~^ + n~ . After some algebra, we can obtain the following 
WPBGK equations for the Wigner components 

^ + ^ Sin kl A-f + 6 . A-/ - er = QVl (3.11) 

H- + ^ sinkl A-/ + bA-f + Cuxf- Qf = Q[f], (3.12) 

e^^ = y(n+ + n--iV^), (3.13) 

whose right hand sides contain collision terms to be described later. Here 

iA^u){x, k) = u{x + 1/2, k) ± u{x - //2, k), (3.14) 

uj = Qq + cDi, (3.15) 

cJo = y (0,0,1), (3.16) 

Q^ = l(0,/3sinfc/A+,27-7cosA;/A+), (3.17) 
n 

b = -(0,f3coskl,'y sinkl), (3.18) 

h 



ejl 

j=-oo 



Qf\x,k,t)= J2 %{F{x,t)),e^^'i;{x,t). (3.19) 



Our collision model contains two terms: a BGK term which tries to send the miniband 
Wigner function to its local equilibrium and a scattering term from the miniband with 
higher energy to the lowest miniband: 

QV] = -—f-^ (3.20) 

n' = ^^±^, n = ^^^^rJ, (3.22) 

2 2 ^ ^ 

^ ik- n±) = — -f- / ^, , ;^ U..M4 \^\l + e^BT dE, (3.23) 



n 



TTh' J_^T^ + [E-SHk)]^ , ^ 

— / (j)^{k;n^)dk. (3.24) 

The chemical potentials of the minibands, /x"*" and /i~ are calculated in terms of n"*" and n~ 
respectively, by inserting (I3.23P in (13.241) and solving the resulting equations. Our collision 
model should enforce charge continuity. To check this, we first calculate the time derivative 
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of n± using (IXTOj) to (l3T2|) : 

^ + ^/ sinkl{f±u.f)dk+'-^ {h.f±V.hf)dk (3.25) 

/A- z"^/' - /A- r/' 

u-ujxfdk^ / u-Qfdk 



n^ 



271" J^-n/l 27r J_^/, 

/A- r^/^ 

= ^^ mn^^-Q[f])dk = T—. 

where we have employed j Qf^dk = 0. Then we obtain: 



— (n+ + n ) + A" 



- / (-smklf + b-f)dk 

T^J-n/l^rl J 



0. (3.26) 



'-tt// 

Noting that A~u{x) = ld{u{x))i/dx, we see that this equation corresponds to charge con- 
tinuity. Differentiating in time the Poisson equation 03. 131) . using 03.261) in the resuh and 
integrating with respect to x, we get the following nonlocal Ampere's law for the balance of 
current: 

dF 



'^ 



(- r (^^inklf + h-pjdk\ =J{t). (3.27) 



Here the space independent function J{t) is the total current density. Since the Wigner 

components are real, we can rewrite 03.27P in the following equivalent form: 

r)F 1p 
^^-y («Im/°-/5ReA2 + 7lmA3>i = J(t). (3.28) 

IV. DERIVATION OF BALANCE EQUATIONS BY THE CHAPMAN-ENSKOG 
METHOD 

In this Section, we shall derive the reduced balance equations for our two-miniband SL 
using the Chapman-Enskog method. First of all, we should decide the order of magnitude 
of the terms in the WPBGK equations 03.111) and 03.12p in the hyperbolic limit. Recall that 
in this limit, the collision frequency 1/r and the Bloch frequency eFMl/h are of the same 
order, about 10 THz for the SL of Section [ITl Typically, 2g/h is of the same order, so that 
the term containing ujq should also balance the BGK collision term. What about the other 
terms? 

The scattering time Tsc is much longer than the collision time r, and we shall consider 
t/tsc = 0(A) <^ 1. Moreover, the gap energy is typically much larger than the miniband 
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widths or the spin-orbit coefficient and a rich dominant balance is obtained by assuming 
that j3/g and 'j/g are of order A. Then we can expand the unit vector z7 as follows: 



z7 = (0, 0, 1) + — sin kl (0, 1, 0) - A^ 
9 



-i sinkin-coskl) (0,1,0) 



(4.1) 



/^^sin^fc/ , 
+ . . (0,0,1) 



'29' 



0(A3 



In this expansion, we have inserted the book-keeping parameter A which is set equal to 1 at 
the end of our calculations (cf. Section [TTl)- From ( 13.1ip and ( 13.12p . we can write the scaled 
WPBGK equations as follows: 



L/-fi = -A [r^ + M 



(4.2) 



Here the operators L and A are defined by 

/ \ 



i^f = f-Tef + 5^ 



\ / 



(4.3) 



Af = 52 







where 



ar I Tb ■ f 

— sinklA-f + A- 



Oi = -^, 62 = — . 







(4.4) 



(4.5) 



The expansion of u in powers of A gives rise to a similar expansion of Q and A. 

To derive the reduced balance equations, we use the following Chapman-Enskog ansatz: 



m=l 



8F °° 

m=0 

00 



dt 



dt 



(4.6) 
(4.7) 
(4.8) 



m=0 



The functions A'^ and Jm are related through the Poisson equation (13.131) . so that 

/ dJm. 



A+ _|_ A- 



e dx 



(4.9) 
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Inserting (gSD to (gSD into (g^), we get 



L/(2) = ^2 - r 



5/(0) 



9t 
9/(1) 



9t 






5/(0) 



9t 



AJW, 



(4.10) 
(4.11) 

(4.12) 



and so on. The subscripts and 1 in the right hand side of these equations mean that we 
replace edF/dt\m = JSom — Jm, dn^/dt\m = A^. Moreover, inserting (14. ip and (14. 6 p into 
(I3.10p yields the following compatibility conditions: 



fr=o, tf^^=^im/r^ 



9 



(2)0 



0, 



(4.13) 
(4.14) 



(2)3 _ /3 ^ ^(1)2 , (3'^ / .(0) 



(0)3n /57. I A0)2 f: 



(0)2^ 



- Im/r^^ + f^ (/o^^^^ - Re/ri - ^ Im f^' - 



9 



^9 



etc. 



To solve (I4.10p for /(") = yj, we first note that 

oo 



]=-oo 



■^,-^w. 



h 



(4.15) 
(4.16) 



Then (HlOD . fl3:22|) and dH]) yield 



^; = ^ 






2 ' 



(4.17) 



where we have used that the Fourier coefficients 



'-^ — ' cos{jkl) (jy^ dk, 



TT 



(4.18) 



are real because 0^ are even functions of k. Similarly, the solution of (14. lip is f^^^ = ip with 

1 - Mi 



^T = -" 






1 + ^] 



(m = 0,3), 

2 



(1 + Z^,-)' + -^l' ' 

2 _ (l + it9j)r|-(5ir] 



(4.19) 
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Here r is the right hand side of (14. lip . The balance equations can be found in two ways. 
We can calculate A'^ for m = 0, 1 by using the compatibility conditions (14.131) and (I4.14p 
in Equations (14.111) and (14.121) . respectively. More simply, we can insert the solutions (I4.17P 
and (I4.19P in the balance equations (13.250 and in the Ampere's law (I3.27p . The result is: 

— - + A-D±(n+,n-,F) = ±i?(n+,n-, F), (4.20) 



dt h 

+j [/3Re(^2^i -alm(^?)i -tMV-?)!] = J, 



n 



"^'^^ Im(V^?±V^3; 



[i + ^i 






+ ^ReV'?±^^^^^^^, (4.22) 



Appendix |X] justifies this second and more direct method by showing that equivalent expres- 
sions are obtained from the compatibility conditions. Note that Eq. (14.210 can be obtained 
from (I4.20p and the Poisson equation. 



V. SPINTRONICS: QUANTUM DRIFT-DIFFUSION EQUATIONS FOR A LAT- 
ERAL SUPERLATTICE WITH RASHBA SPIN-ORBIT INTERACTION 

In the simpler case of a LSL with the precession term of Rashba spin-orbit interaction 
(but no intersubband coupling), we can obtain explicit rate equations for n^ by means of 
the Chapman- Enskog method. In the Hamiltonian (13.10 . we have 7 = (? = 0, so that h^ = 
and z7 = (0, 1,0). However, the Fermi-Dirac distribution is different from (12.60 for a LSL. 
We have to replace En instead of h'^k^/(2m*), sum over n for all populated QW energy 
levels and integrate over ky only. Provided only Ei is populated, we obtain the following 
expression instead of (I3.23p : 

^-00 1 + exp ' ^^^ 



ksT 

where the broadened density of states is 

DAE) = J- r ik,ir m-E\=^r ,e, ^£%=^. (5.2) 
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Note that (E^D becomes the ID density of states D{E) = V2m*e{E)/{2TThL^y/E) as T ^ 
0+. We have not included a factor 2 in (15.21) because all the electrons in each of the 
minibands (with energies S^{k)) have the same spin. Inserting (12.81) in (15.21) and using the 
residue theorem to evaluate the integral, we obtain 



DriE) 



AirhL, 



(5.3) 



X 



+ 



^/WT72^E + f^ + E + ^ - J y^WW^TE + f^ - E - ^, 



V2 



V2 



Ve^ + V2TE + r2 



Ve'^ - V2TE + T^ + E 



^ + VV^^-v^ri5; + P-i^ + ^ 



v/^E2^V2fETr2 



As E ^ +00, Dr{E) ~ V2^/{27rhL,y/E), whereas Dr{E) = 0{\E\-^/^) as E ^ -oo. 
Then the convolution integral (15. ip is convergent. 

In the present case, minibands correspond to electrons with spin up or down which have 
different energy. Scattering between minibands is the same as in (I3.2ip . —{i^f^ + f)/Tsc which 
yields dn^/dt+ . . . = ^n^ /ts^ in (13.251) . only if the chemical potential of the miniband with 
lowest energy, /i~, is less than the minimum energy of the other miniband, E^^^ =in.m].£^ {k) . 
Otherwise (yU~ > ^^in)' ^'^^ scattering term should be — 2//rsc, which yields dri^ /dt + . . . = 
^{n^ — n~)/Tsc in (13.250 . thereby trying to equahze n"*" and n~; cf. Ref.— . 

Now we shall derive the balance equations in the hyperbolic limit using the Chapman- 
Enskog method as in Section [IVl In the scaled WPBGK equations (14. 2p . the operators L 
and A are 




\ aT 

+ -^smklA~f 

2/ + (^/°-/)^(Cin-/^-^ ' ^ 



+ — cos kl A 
n 







+ — sin kl A^ 

h 



( \ 

P 



V-/V 



(5.4) 
(5.5) 



where ^2 is given by (14.51) . d{x) is the Heaviside unit step function and VP = (0+ + )/2, 
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n = (0, 1, 0) ((/)+ - (j)')/2. The hierarchy of equations (HlHl) - KW\ is simply 






5/(0) 



at 

5/(1) 



A/ 



(0) 



at 



A/ 



(1) 



5/(0) 



at 



(5.6) 

(5.7) 

(5.8) 



and so on. The compatibihty and solvabihty conditions are: 



/, 



(m)0 _ f-im)2 







Jo 







(L/('")°)o = (L/(™)- 



0, m > 1. 



(5.9) 



The solution f^^^ = Lp of (15.61) is 

n 0+ + 07 1 - Z^, 



<^i 



^TT#' ^-^-'^ ^?-^^ 



1 - it9j- 



(5.10) 



where we have used that the Fourier coefficients • are real because 0^ are even functions 
of k. Similarly, the solution of (15.71) is /'•^-' = -0 with 

1 - i'd, 



c = ^r 



3 'J 1+^2 



(m = 0,2), V.— V'.— O. 



(5.11) 



Here r is the right hand side of (15.71) . The balance equations can be found in two ways. We 
can calculate A^ for m = 0, 1 by using the solvability conditions (15.91) in Equations (15. 7p 
and (15. Sp . respectively. More simply, we can insert the solutions (I5.10p and (15. lip in the 
balance equations (13.251) and in the Ampere's law (I3.27p . In both cases, the result is: 



at 

dF 

e— + e{D+ + D_), = J, 

D± = -^ A-Im(y,? ±y^l + ij',± ij!) ± ^ A-Re(^? ± V'? + ^? ± i^f) 
n n 



R 



n+ - n- 9{f^- - £+ J 



(5.12) 
(5.13) 
(5.14) 
(5.15) 



NONLINEAR ELECTRON AND SPIN TRANSPORT IN SUPERLATTICES 



21 



A straightforward calculation of fl5.14p yields 

(g^i ± P)^f r (0+ - 0r) [2c^^i ±Pil- ^i)] 

[2a^i ± p{l - dl)]aT d^^ \ ^_ f a^i±P 



(5.16) 



^2(1+ ^2)2 5^± 



A- 



/i(l + ^?) 



h 



a (3^2 _ 1) ± ^^^(3 _ ^2) i^2^± ^j 

^1(1 + ^2)3 



1/ 1 



he \e 
2^2(1 + ?92) 



± (n^ 

CtTsc 

a(0+ + 0r)^i 



n 



1/ 1 



A"n^ 



r 



2^2(1 +^2^ 



(a 



P^T^aP^i)^' 



^f 



1+^2 



+ 



a^)^! T 2a/3] A" 



^20f 



1 



^Lj 



We have numerically solved the system of equations (15.121) - (15.161) . with the following 
boundary conditions in the interval —21 < x < 0: 

dF 



'm^^^ 



J, 



n 



n 



Nn 
2 '' 



whereas in the collector Nl < x < N {I + 2), ( I5.17P and 



dx 







(5.17) 
(5.18) 

(5.19) 



hold. We have used the following values of the parameters: a = Ai/2 = 8 meV, /3 = 2.63 
meV, dw = 3.1 nm, ds = 1.96 nm, / = dw + ds = 5.06 nm, L^ = 3.1 nm, T = 5 K, 
T = 5.56 X 10-1^ s, Tse = 5.56 x IQ-^^ g^ ^^ ^ 4948 x 10^° cm-2, m* = {0m7dw + 
0.15dB)mQ/l, V = 3V, N = 110. We have used a large conductivity of the injecting contact 
a = 11.78 fi^^m^^. With these values, we select the following units to present graphically 
our results: Fm = h/^elr) = 23.417 kV/cm, xq = eFMl/ieND) = 19.4 nm, to = h/a = 0.082 
ps, Jo = aeND/i2h) = 3.94 x 10^ A/cm2. 

Fig. [2](b) - (d) illustrates the resulting stable self-sustained current oscillations. They 
are due to the periodic formation of a pulse of the electric field at the cathode x = and 
its motion through the LSL. Fig. MJo) depicts the pulse when it is far from the contacts 
and the corresponding spin polarization is shown in Fig. M^d). It is interesting to consider 
the influence of the broadening F and the Fermi-Dirac statistics on the oscillations. At 
high temperatures, Boltzmann statistics and a semiclassical approximation should provide 
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a good description. The semiclassical approximation is equivalent to dropping all spatial 
averages in our previous formulas. Since Xq ^ I, the effect of dropping spatial averages 
should be rather small. Using Boltzmann statistics yields explicit formulas for /x^ in terms 
of n^. In fact, we only have to replace e^^ -E)/{kBT) j^g^g^d of the 3D Fermi distribution 
[1 + e^^~^ )/('=sT')j-i jj^ j7;q_ fl5.ip . Using the relation fl3.24p between n^ and 0^, we obtain 



/ a cos fci=p/3 I sinfcZI 

n± ''"''^^ , '"^ ^, (5.20) 

j;dKeicY> r"°"TBT 



and therefore, 



b^ = n^ -, ^ ^ '-, (5.21) 



lo dK exp ^^ ^^y 

for J = 0, 1, . . . Similar relations hold for the case of a SL with Boltzmann statistics in the 
tight-binding approximation. 

The results are shown in Fig. [2j Fig. [2]^a) depicts the relation between electron current 
and field for a spatially uniform stationary solution with n^ = Nd/2. We observe that all 
curves are similar. However the curves for F = and F = 1 meV are close while the curve for 
F = 5 meV has dropped noticeably. The shapes of J{t) for F = and F = 1 meV in Fig. [2t^b) 
are close and quite different from that for F = 5 meV. If we look at the corresponding field 
profiles in Fig. [2](c) and (d), for F = and F = 1 meV the oscillations of the current are 
caused by the periodic nucleation of a pulse of the electric field at x = and its motion 
towards the end of the LSL. The pulse far from the contacts shown in Fig. [2](c) is larger in 
the case of F = than for F = 1 meV. In the case of F = 5 meV (not shown), the pulse 
created at x = becomes attenuated and it disappears before arriving at x = Nl. This 
seems to indicate that the lowest voltage at which there exist stable self-sustained current 
oscillations is an increasing function of F: If we fix the voltage at 3 V and increase F, the 
critical voltage threshold to have stable oscillations approaches our fixed voltage of 3 V. 
Then the observed oscillations are smaller and the field profiles correspond to waves that 
vanish before reaching the end of the device, as it also occurs in models of the Gunn effect 
in bulk semiconductors^. 
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S^^ 



(c) 



r = o 

r = 1 meV 
r = 5 meV 
Boltzmann statistics 




14 16 18 



F/F 



M 



20 22 24 26 28 




1500 2000 




FIG. 2: (a) Electron current vs field in a spatially uniform stationary state for different values of 
the broadening T using the Fernii-Dirac distribution and for the Boltzmann distribution without 
broadening, (b) Total current density vs time, and the (c) electric field and (d) spin polarization 
profiles during current self-oscillations for F = (solid line) and 1 meV (dashed line). Parameter 
values are N = 110, No = 4.048 x 10^° cm"^ ds = 1.96 nm, L^, = dw = 3.1 nm, I = 5.06 nm, 
r = 0.0556 ps, t^c = 0.556 ps, F = 3 V, cr = 11.78 f^-^m-i T = 5 K,ma = 8 meV, /? = 2.63 meV. 
With these values, Ai = 16 meV, xq = 19.4 nm, to = 0.082 ps, Jq = 3.94 x lO'' A/cm^. 

VI. CONCLUSIONS 



We have presented a Wigner-Poisson-BGK system of equations with a collision broad- 
ened local Fermi-Dirac distribution for strongly coupled SLs having only one populated 
miniband. In the hyperbolic limit in which the collision and Bloch frequencies are of the 
same order and dominate all other frequencies, the Chapman-Enskog perturbation method 
yields a quantum drift-diffusion equation for the field. Numerical solutions of this equation 
exhibit self-sustained oscillations of the current due to recycling and motion of charge dipole 
domains^. 

For strongly coupled SLs having two populated minibands, we have introduced a peri- 
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odic version of the Kane Hamiltonian and derived the corresponding Wigner-Poisson-BGK 
system of equations. The colhsion model comprises two terms, a BGK term trying to bring 
the Wigner matrix closer to a broadened Fermi-Dirac local equilibrium at each miniband, 
and a scattering term that brings down electrons from the upper to the lower miniband. By 
using the Chapman-Enskog method, we have derived quantum drift-diffusion equations for 
the miniband populations which contain generation-recombination terms. As it should be, 
the recombination terms vanish if there is no inter-miniband scattering and the off-diagonal 
terms in the Hamiltonian are zero. These terms may represent a Rashba spin-orbit interac- 
tion for a lateral superlattice. For a lateral superlattice under dc voltage bias in the growth 
direction, numerical solutions of the corresponding quantum drift- diffusion equations show 
self- sustained current oscillations due to periodic recycling and motion of electric field pulses. 
The periodic changes of the spin polarization and spin polarized current indicate that this 
system acts as a spin oscillator. 
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APPENDIX A: BALANCE EQUATIONS FROM COMPATIBILITY CONDI- 
TIONS 

We know that ip^ = ip'^ = from (14. lip . Then the compatibility conditions fl4.13p and 
fl4.14p become 



^° = 0, ^o' = 0, 






(A.l) 


fr=o, fr- 


= ^-^^'+4^^^^°- 


-Rev^D, 


(A.2) 



Equations (lA.ip imply that (L^)" = for m = 0, 3 in (HTTl) . Since ip^ = (n+ + n")/2 and 
ifl = {n^ — n^)/2, these conditions yield 



2 dt 



-- A"Imy?? - ^ A-lmcpl = 0, 
n n 



2 dt 



( y i V / 
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wherefrom we obtain 

A^ = T- + ^A-lm{^1±^l). (A.3) 

Tsc "- 

Let us now calculate Af. Equations ([Ql) imply (L/(2))0 = and (L/(2))3 = f^"^^^ given 
by ([QD in KW\ . After a little algebra, we find 

Af = ^^4^ A-Im(^° ± i^l) - ^ (A-Re^i^ ^ A+Im^J) (A.4) 



gr 8g^T 



T^ limPl ± ^ [2Re<^^ + 0+ - ^s" - 2(n+ - n-)]. 



We will now transform (1A.4P in an equivalent form by eliminating Iieip2 and Im^^ in 
favor of Reipl and Im^^, respectively. Eq. (14.101) implies that (1 + i^2)'^2 = i4>2 ~ ^2)/'^^ 
and therefore, 

Re^l = ^2lm^l + ^l^^. (A.5) 

Similarly, Eq. (14.111) implies that (1 + i'di) ipf + Siipl = rf, and therefore, 

ImV'i = -t9i Re^pf - 5i Im^l + Imrj. (A. 6) 



The right hand side of 


(41U) yields 






















rl- 


2<7 


n-e--- 


' — 


0" 




/5r 


A 

V 


14 


. ^-i2kl 
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^ 

J 






\ ll 





wherefrom 






























Imr? - f (0^ 


02 


" — 


n+ + 


n-) 


/5r 
2n 


A~ 


Imyjg- 







(A.7) 



Inserting (lA.Sp . (]A.6P and (]A.7p in (]A.4p . we obtain the equivalent form: 



Af = ^^^ A-Im(^° ^ ^^i) ~ f iA-Re^Pl ± A+Imz/'l) (A.8) 

2/3 /3 /3^ /3^ 

±4- ImV'l ± — ?9iRe^? ± — ;- ^2 Im<^2 ± 4- A-Im(^°. 
n (jfT 4(7"'r 2ng 

Inserting ([Qj) and this expression in (jM]) and using (KIT} , yield (g^Q]), (l4:22|) and (14:231) . 



Up to order A^, we have thus proven the following statement: 

By using the compatibility conditions in the hierarchy of equations l[4-ll }, (4jJ^> "^^ 



obtain the same balance equations for n as by direct substitution of the solutions of the 
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hierarchy into equations Ii3.25\) (which arise from integration of the kinetic equation over k). 
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